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Factorization of Difference Equations by Semiconjugacy 
with Application to Non-autonomous Linear Equations 



H. SEDAGHAT* 



Abstract. The existence of a semiconjugate relation permits the transformation of a higher 

f^ ' order difference equation on a group into an equivalent triangular system of two difference 

Cn , equations of lower orders. Introducing time-dependent form symmetries in this paper enables 

^~i| us to identify the semiconjugate property in a larger set of non- autonomous difference equations 

^^ ' than previously considered. We show that there is a substantial class of equations having this 

feature that includes the general (non-autonomous, non-homogeneous) linear equation with 
C^ , variable coefhcients in an arbitrary algebraic field. 

, , 1 Introduction 

P^ . Difference equations of order greater than one that are of the following type 

•i-H ■ 

determine the forward evolution of a variable x„ in discrete time since the time index or the 
^ ' independent variable n is integer-valued with n > 0. 

^ . In previous studies of semiconjugate factorizations of difference equations of type ([I]), e.g., [3], 

tH- ! |1]) [5] or [6], the form symmetry linking the higher dimensional unfolding map of the original 

CN ' equation to that of the lower dimensional factor was assumed to be independent of n. While 

l/^ . this assumption did not substantially curtail the applicability of the method, it did rule out certain 

^^ I non-autonomous equations. For example, the method worked for non- homogeneous linear equations 

with constant coefficients but did not apply to linear equations with variable coefficients. 

The main goal of this article is to extend the aforementioned factorization method to allow time- 
dependent form symmetries where the form symmetry may depend explicitly on the independent 
K> I variable n. This extension is significant as it covers all non- autonomous equations of type ([T]). In 

C^ ■ particular, the extended method may be applied to general (non- autonomous, non-homogeneous) 

linear equations over arbitrary algebraic fields to show that such equations admit semiconjugate 
factorizations via eigensequences (i.e., the solutions of an associated discrete Riccati difference 
equation of lower order). For ease of reference we state some of the basic concepts and notation 
here; additional background material for this article is available in [5J. 
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As usual, the number A: in ([T]) is a fixed positive integer and k + 1 represents the order of the 
difference equation ([TJ. The underlying space of variables Xn is a group G and /„ : G'^+^ — )■ G 
is a given function for each n > 1. If /„ = / does not expHcitly depend on n then ([T]) is said 
to be autonomous] it is non- autonomous otherwise. A (forward) solution of Eq.([T]) is a sequence 
{xn}'^^_k that is recursively generated by ([T]) from a set of A; + 1 initial values xq, x_i, . . . , x„fc G G. 
Forward solutions have traditionally been of greater interest in discrete models that are based on 
Eq.([T]) although other types of solutions (e.g., those having domain Z, the set of all integers) can 
also be readily defined. 

Each /„ is "unfolded" by the associated vector map F^ : G^^^ — t- G^^^ that are defined as 

Fn{uo,...,Uk) = [fniuo,... ,Uk),Uo,... ,tifc-l], Uj £ G ioi j = 0,1, . . . , k. (2) 

The unfoldings F^ determine the equation 

(yO,n+l, yi,n+l, • • • , yk,n+l) = FniUQ^n, 2/1,71, • • • , yk,n) 

in G'''^^. Each vector (j/o.n+i, ■ ■ ■ ■, yk,n+i) represents a state of the system, or of Eq.([l]); G^^^ is the 
state space, in analogy to the phase space in differential equations. 

2 Semiconjugate relation and factorization 

Let Fn be the unfolding on G'^^^ of /„ for each n. Then ([1]) is equivalent to 

Xn+l=Fn{Xn), Xn = {Xn, ■ ■ ■ , Xn-k)- (3) 

We are interested in deriving a lower dimensional equation 

y„+l = ^niXn), Yn = (Vn, • • • , Vn-m+l), m < k (4) 

for ([3|). If there exists a sequence of maps Hn '■ G^^^ — )• G*" such that for every solution {X„} of 

m 

Yn = Hn{Xn), n = 0,1,2,... (5) 

is a solution of (jl]) then 

^n{Hn{Xn)) = ^n{Yn) = Y^+l = F„+l(X„+i) = /7„+i(F„(X„)). 

Therefore, ^ is satisfied for all solutions of ([3]) and (JH if and only if the sequence {Hn} of 
maps satisfies the following equality for all n 

Hn+l oFn = ^nOHn. (6) 

If the mappings Hn are independent of n, i.e., Hn = H for all n then Eq.([6]) reduces to 

HoFn = ^nOH (7) 

namely, the time-independent semiconjugate relation as defined in prior studies. We can now give 
the following more general defintion. 



Definition 1 Let k > 1, 1 < m < k. If there is a sequence of surjective maps Hn : G^^ — )• G^ 
such that Eq.(^ is satisfied for a given pair of function sequences {Fn} and {^n} then we say that 
Fn is semiconjugate to ^n for each n and refer to the sequence {Hn} as a (time- dependent) 
form symmetry of Eq.^ or equivalently, of Eq.^^. Since m < k + 1, the form symmetry {Hn} 
is order-reducing. 

Technically, a time-dependent form symmetry can also be defined as a single map 

H-.Nx G^+^ -^ G"*, H{n- uq, . . . , Uk) = i^„(uo, . . . , n^) 

We choose the sequence definition due to its more intuitive content. 

The following result extends its time-independent analog in [5] and makes precise the concept 
of semiconjugate factorization for the recursive difference equation ([1]). 

Lemma 2 Let k > 1, 1 < m < k, let h^ ■ G ''^'^ — )■ G for n > —m + 1 be a sequence of functions 
on a given non-trivial group G and define the functions Hn : G^^ — )• G™" by 

Hniuo,...,Uk) = [no * hn{ui, . . . ,Uk+l-m), ■ ■ ■ ,Um-l * K-m+liUm, ■ ■ ■ ,Uk)]. (8) 

Then the following statements are true: 

(a) The function Hn defined by ^ is surjective for each fixed n > 0. 

(b) If {Hn} is an order-reducing form symmetry then the difference equation ([7P is equivalent 
to the system of equations 



tn+l — ^^n\J"ni ■ ■ ■ i ^n— m+1 j; V>^j 

a^n+l = tn+l * hn+l{Xn, ■ ■ ■ , Xn-k+m)~ (10) 

whose orders m and k -\-\ — m respectively, add up to the order of ([I]j. 

(c) The map <!>„ : G™' — >■ G™' in ^B^ is the unfolding of Eq.l§^ for each n > 0; i.e., each $„ is 
of scalar type. 

Proof, (a) Let n be a fixed non-negative integer and for j = 0, . . . , ?7i — 1 denote the j'-th coordinate 
function of Hn by 

r]j+l{uo,...,Uk) =Uj * hn-jiuj+i, . . . ,Uj+k+l~m) (11) 

Now choose an arbitrary point {vi, . . . , Vm) € G™" and define 

Um-l =Vm* hn-m+l{Um,Um+l •••,Mfe)~"^, (12) 

Um = Um+1 = ...Uk=U 



where u is a fixed element of G, e.g., the identity. Then 

Vm = Um-l * hn-m+l {u,U...,u) 

^ ^m— 1 * "n— m+1 ["Umi ^m+1 • • • ) "^fc j 

= Vmiuo,- ■ ■ ,Uk) 

= 1]m{uO,- ■ ■ ,Um-2,Vm * hn-m+l{u,U . . . ,u)~'^,U. . . ,u). 
^ . ' 

Um-l 

for any selection of elements uq, . . . , Um-2 S G. Using the same idea, define 

Um-2 = Vm-1 * hn-m+2{Um-l,U . . . , u)~ 

with Um-i defined by ([12]) so as to get 

Vm-l = Um~2 * K-m+2iUm-l,U . . . ,u) 

= Um-2 * hn-m+2{Um-l,Um ■ ■ ■ ,Mfc-l) 

= rjm-iiuo,...,Uk) 

= 11m-l{uo,- ■ ■ ,Um-3,Vm-l * hn-m+2{Um-l,U . . . , 'u)~\ ^^-1 , « • • • ,u) 
^ V ' 

Um-2 

for any choice of uq, . . . ,Um-3 G G. Continuing in this way, by induction we obtain elements 
Um-i, ■ ■ ■ ,uo G G such that 

Vi = r]i{uo,...,Um-i,u...,u), i=l,...,m. 

Therefore, Hn{uo, . . . , Um-i,u . . . ,u) = {vi, . . . , Vm) and it follows that Hn is onto G*". 

(b) To show that the SC factorization system consisting of equations Q and p^ is equivalent 
to Eq.([T]) we show that: (i) each solution {xn} of ([I]) uniquely generates a solution of the system ([9]) 
and ([To]) and conversely (ii) each solution {{tn, Vn)} of the system ([9]) and ([TO]) correseponds uniquely 
to a solution {x„} of ([T]). To establish (i) let {x„} be the unique solution of ([T]) corresponding to 
a given set of initial values xq, . . . x_i^ G G. Define the sequence 

t-n ^ Xn * iln\Xn—l) • • • ) Xn—k+m—l) \^'^) 

for n > —m + 1. Then for each n > if Hn is defined by ([8]) it follows from the semiconjugate 
relation ([6]) that 

Xn+l — Jn\Xni • • • ) Xn—k) 

— (pny^n * i^n\Xn—lj • • • ; ^n— fc+m— 1 j) ■ ■ ■ j 

Xn—m+1 * i^n—m+\ \Xn—rm ■ ■ ■ i ^n—kj) * ['^n+1 y^m • • • ; Xn—k+m)\ 

— Vny^nt • • • i t-n— m+1 j * [P'n+lyXni ■ ■ ■ i Xn—k+m)\ 



Therefore, (/»„(*„,..., tn-m+i) = Xn+i * K+iixn, • • • , Xn-k+m) = tn+1 so that {tn} is the unique 
solution of the factor equation ([9|) with initial values 

t—j = X—j * /i_j(x_j_x, . . . , a;_j_fc_|_m_xj, j = (J, . . . , ?7i — 1. 

Further, since Xn+i = tn+i * [hn+i{xn, ■ ■ ■ ,Xn-k+m)]^^ for n > by (fT3|). {x„} is the unique 
solution of the cofactor equation (|1U|) with initial values y_j = x^j for i = 0,1, . . . ,k — m and with 
the values t„ obtained above. 

To establish (ii) let {(t„,yn)} be a solution of the factor-cofactor system with initial values 

*0) • • • ,t-m+l,y-m, ■ ■ -V-k € G. 

Note that these numbers determine y_m+i, ■ ■ • ,2/o through the cofactor equation 

y-i = *-i * [^-i(y-i-l'•••'y-i-l-fc+n^)]"^ j = 0, ...,m-l. (14) 

Now for n > we obtain 

Vn+l = tn+1 * [hn+liUn, ■ ■ ■ , ?/n-fe+m)] 

^ (Pn\tm • • • ) tn—m+l) * ["n+liym • • • ) yn~k+'m)\ 
= (t>n{yn * hniVn-l-,- ■ ■ , 2/n-A:+m-l), • • • , 

yn—m+l * iT'n—Tn+iyyn—mi • • • ) yn~k)) * "ra+liyn; • • • i yn—k+m) 
^ Jn\ym • • • ) yn—k) 

Thus {2/n} is the unique solution of Eq.([T]) that is generated by the initial values (fHl) and 
y_m,, • • • y-k- This completes the proof of (b). 

(c) We show that each coordinate function (j)j^n is the projection into coordinate j — 1 for j > 1. 
From the definition of Hn in ([8]) and the semiconjugate relation ([6]) we infer that 

Hn+l{Fn{uo, . . . ,Uk)) = Hn+l{fn{uo, ■ ■ ■ ,Uk),Uo,. ■ ■ ,Uk-l) 

= {fn{uO,--- ,Uk) * K+l{uo,... ,Uk-m), 
Uo * hn{ui,... ,Uk-m+l),--- , 
Um-2 * hn-m+2{Um-l, ■ ■ ■ ,Uk-l))- 

Matching the corresponding component functions in the above equality for j > 2 yields 

4>j,n{uo * hn{ui,.. . ,nfc+i_m),... ,Um-l * /ln-m+1 (-"m, • • • ,Uk)) = 
Uj-2 * h{uj^i,Uj ... , Uj+k-m-l) 

which shows that (pj^n maps its j-th coordinate to its (j — l)-st. Therefore, for each n and every 
{ti,...,t^)e Hr,{G^+^) we have 

^nih,- ■ ■ ,trn) = [(pnih, ■ ■ ■ ,tm),tl, ■ ■ ■ ,tm-i] 



i.e., ^n|H„(G*+i) is of scalar type. Since by Part (a) Hn{G^^^) = G^ for every n, it follows that 
$„ is of scalar type. ■ 

The pair of equations Q and pUj) in Theorem[2]is uncoupled in the sense that ([S]) is independent 
of (jlOp . Such a pair forms a triangular system as defined in [1] and [7]. In the next definition we 
use convenient and suggestive terminology to describe these equations. 

Definition 3 Eq.^ is a factor of Eq.^^ since it is derived from the semiconjugate factor <!>„. 
Eq.UO\) that links the factor to the original equation is a cof actor of Eq.(l^. We refer to the 
system of equations ^ and ilO\) as a semiconjugate (SC) factorization of Eq.m). Note that 
orders m and k + 1 — m of (0j and < [J0|) respectively, add up to the order of |7p. We refer to the 
system of equations (0j and dJOj) as a type-(m, k + 1 — m) order reduction of Eq.m). 

3 Invertible-map criterion 

In [4] and [5] a useful necessary and sufficient condition is obtained by which to determine whether 
the difference equation ([T]) has order-reducing form symmetries (not time-dependent). In this 
section we show that the same useful idea extends to the time-dependent context. Applications 
and examples are discussed in the next section. 

Consider the following special case of ([8]) with m = k 

Hn{uo,Ui,... ,Uk) = [uq * hn{ui),Ui * /l„_i (^2), . . . , "Ufc-l * K-k+l{uk)] (15) 

with hn : G ^ G being a sequence of surjective self-maps of the underlying group G for n > 
—k + 1. If (dJ has the form symmetry (jlSp then it admits a type-(A;, 1) order-reduction and its SC 
factorization is 

tn+l = (f'nitn, • • • , ^n-fc+l), (16) 

Xn+1 = tn+1 * hn+l{Xn)~^ . (17) 



The initial values of the factor equation (I16p are 

t-j = x-j * h_j{x-j+i), j = 0, 1, . . . , fc - 1 

Theorem 4 (Time- dependent invertible map criterion) Assume that hn ■ G ^ G is a sequence of 
bijections of G for n > —k + 1. For arbitrary elements uo,vi, . . . ,Vk G G and every n > define 
Co,n(no) = uo and for j = 1,..., k, 

Cj,niuO,Vl,...,Vj) = h:^lj^-^{Cj-l,n{uo,Vl, . . . ,Vj-iy'^ * Vj). (18) 

with the usual distinction observed between map inversion and group inversion. Then Eq.^ has 
the form symmetry {Hn} defined by i f75|) if and only if the quantity 

/n(Co,n, Cl,ni^O,Vl), . . . , Ck,niuo, Vi, . . . , Vk)) * hn+liuo) (19) 



is independent of uq for every n > 0. 

In this case Eq. ^ has a SC factorization whose factor functions in / tJ6]) are given by 

4>n{vi,- ■ ■ ,Vk) = fn{Co,n,Cl,niuo,Vl),... , Ck,n{uo,Vl, . . . ,Vk)) * hn+l{uo) . (20) 

Proof. Assume first that (J19p is independent of uq for all wi, . . . , f^ so that the functions 

(l)nivi,...,Vk) = fniCo,n,Cl,n,---,Ck,n) */ln+l(?io) (21) 

are well defined. Next, if Hn is given by (fT5|) then for all uq, ui, . . . , Ufc 

(l>n{Hniuo,Ui,... ,Uk)) = (/'„(uo * hn{ui),Ui * /l„_i (-0,2), • • • , "U-fc-l * K^k+liuk)) ■ 

Now, by (jlSp for each n and all mQ) "Wi 

Cl,n{'^0,UO*hn{ui)) = K^ {Uq^ * Uq * hn{ui)) = Ui- 

Similarly, for each n and all uo,ui,U2 

C2,ni'^0,Uo *hniui),Ui */l„_i(u2)) = h~^^{Cl,n{uO,Uo * hn{ui))~^* 

Ul *hn-l{u2)) 
= U2. 

Suppose by way of induction that 

C«,n(^0 * hn{ui),... ,Ul^i * hn-l+l{uk)) = Ul 

for 1 < ^ < i- Then 

Cj,n{uo * hniui), . . . ,Uj-i * hn-j+l{Uj)) = k^^ j _^_-^^{uj\ * Uj-i * hn-j+l{Uj)) = Uj . 

Thus by ^ 

4>n{Hn{uo,Ui,... ,Uk)) = /„('Uo, . . . , Ufc) * hn+l{uo) 

Now if Fn and $„ are the unfoldings of fn and (pn respectively, then 

Hn+l{Fn{uo,. . . ,Uk)) = [fn{uo,- ■ ■ ,Uk) * /in+l (uq), ""-O * hn{ui), 

. . . , Uk~2 * hn-k+2{uk~-l)] 
= [4>n{Hn{uo,Ui,... ,Uk)),Uo * hn{ui), 

..., Uk-2 * hn-k+2{Uk-l)] 
= ^n{Hniuo,...,Uk)) 



and it follows that {Hn} is a semiconjugate form symmetry for Eq.([T]). The existence of a SC 
factorization with factor functions defined by (j20p now follows from Lemma [2j 

Conversely, if {H^} as given by psp is a time-dependent form symmetry of Eq.([T]) then the 
semiconjugate relation implies that for arbitrary uq, . . . ,Uk in G there are functions (/>„ such that 

/„(uo,...,Mfc) *hn+l{uo) = (j)niuo * /l„(lii), . . . , Ufc.i * /l„_fc+i(Ufc)). (22) 

For every uo,vi, . . . ,Vk in G and with functions Cj,n as defined above, note that 

Cj-l,niuo,Vi,...,Vj-i) * hn^j+l{Cj,niuO,Vl, . . . ,Vj)) = Vj, j = 1, 2, . . . , /c. 

Therefore, abbreviating Cj,niuo,vi, . . . ,Vj) by Cj,n we have 

fn{Co,n,Cl,n,--- ,Ck,n) * hn+l{uo) = 0n(Co,n *^n(Cl,n),Cl,n * /in-l(C2,n), 

• • • 5 Cfc-l,n * hn-k+l{Ck,n)) 

= 0„(7;i,... ,t;fc) 

which is independent of uq. m 

Recall that an algebraic field T ={T, +, •) is, in particular, a commutative group with respect to 
addition. Further, its set of nonzero elements J-"\{0} is a commutative group under multiplication. 
A simple yet important type of form symmetry may be defined on a field. 

Definition 5 Let T he a non-trivial field and {««} o. sequence of elements of J- such that an £ 
J^\{0} for all n > —k + 1. A (time- dependent) linear form symmetry is defined as the following 
special case of [TB\) with hn{u) = —an-iu 

[Uo - Q„_l-Ui, -Ui - an-2U2, ■■■ , Uk-1 - ttn-kUk]- (23) 

The sequence {an} of nonzero elements in T may he called the eigensequence of the linear form 
symmetry. If Eq.^^ has a linear form symmetry then {an} is an eigensequence of (QP. 

The use of the term "eigen" which is borrowed from the theory of linear equations is apt here 
for two reasons. First, the sequence {an} characterizes the linear form symmetry (j23p completely 
and secondly, we find below that linear difference equations indeed have linear form symmetries. 

The existence of a linear form symmetry implies a type-(A:, 1) order reduction for Eq.([T]) and a 
SC factorization where the cofactor equation (J17p is determined more specifically as 

a^n+l = tn+l + anXn- (24) 

The following necessary and sufficient condition for the existence of a time-dependent linear 
form symmetry is an application of Theorem^ We drop further mention of "type-(A;, 1)" as we do 
not discuss any other order reduction types in the remainder of this paper. 



Corollary 6 Equation |ip has a time- dependent linear form symmetry of type Ii23\) with an eigense- 
quence {an} in a non-trivial field J- if and only if the quantity 

fn{uO,Cl,niuo,Vl),... , Ck,n{uo,Vl, . . . ,Vk)) - ttnUQ (25) 

is independent of uq for all n >0 with the functions Qj^n for j = 1, . . . ,k given by 

Cj-l,niuo,Vi,...,Vj-i) -Vj 



0,n(no,fl,--- ,Vj) 



Ui=i «"- 



3 ^ 

i=l p=l 



Proof. The conclusions follow immediately from Theorem |4] using hn{u) = —an-iu. The last 
equality above is established from the equality preceding it by routine calculation. ■ 

Remark 7 // Eq. (Cp has a linear form symmetry then by Corollary an eigensequence of (OP can 
be defined equivalently as a sequence {an} in J^\{0} for which the quantity in \25\) is independent 
of uq for every n> 0. 

We close this section with an example of a nonlinear equation that has a linear form symmetry. 
For additional results and examples, we refer to [2]. 

Example 8 Consider the following third-order nonlinear difference equation 

Xn+l = (-1)"+-^X„ - 2x„_i + gn{Xn + Xn-2) (26) 

where (7^ : M — )• R zs a given function for each n > —2. By Corollary IS\ a linear form symmetry for 
/i26\) exists if and only if the quantity 



(-l)"+lno - 2Cl,n + 9n{uo + C2,n) " OnUQ (27) 

is independent of uq for all n. Substituting 

^ UO-Vl Cl,n -V2 Uo-Vi- a„__if2 

(,l,n — , C,2,n — — • 

On-l 0!n-2 Oin-ian-2 



in (21) and rearranging terms gives 



(-1)"+' - ^- - an") no + -^vi + 

5n ( ( IH ) ""0 Vi V2 

VV an-ian-2j Oin-\an-2 OLn-2 



which is independent of uq for all n if the coefficients of the uq terms are zeros; i.e., for all n, the 
numbers an satisfy both of the following equations 

an = (-1)"+^ - ^- (28) 

0!n-l 

an-1 = —. (29) 

0!n-2 

Every solution of Eq. (2§\) is a sequence of period 2 

Q, — ,Q, — ,■■■> (30) 



where q = q_2 G M. Now [29^) yields a_i = —1/(7, which we substitute as an initial value in Eq.[2S\) 
to get 

ao = -1 + 2q. 

Now to make the period-two sequence in i30\) also a solution of i28\} . we require the above value of 
ao to be equal to q; thus 

ao = q^2q— l = q^q = l. 

We check that if ao = q = 1 in [2^) then 

2 1 2 

«! = 1 = — 1 = — , 02 = —1 = 1 = q, etc 

q q -1 

so that both of the equations 112^) and 1129]) generate the same sequence {««} where an = (—1)" for 
n > —2. It follows that {(—1)"} is an eigensequence for i26\). 

4 Factorization of linear equations 

We expect that linear difference equations are among difference equations that have the hnear form 
symmetry and this is indeed the case. The following application of Corollary [6] and Theorem U] 
gives the semiconjugate factorization for non-autonomous and non-homogeneous linear difference 
equations. 

Corollary 9 (The general linear equation) Let {ai^n}, i = 1,. ■ ■ ,k and {bn} be given sequences 
in a non-trivial field J- such that a^^n 7^ for all n > 0. The non-homogeneous linear equation of 
order fc + 1 

Xn+l = ao,nXn + ai^nXn-l H h Ok-nXn-k + K (31) 

has a linear form symmetry with eigensequence {a„} for every solution {««} in T of the fol- 
lowing Riccati equation of order k 

an = ao,n -\ \ \ \ (32) 

an-1 Oin-ian-2 "n-l • • • «n-fc 

10 



The corresponding SC factorization of h31\l is 

k k 

tn+l = ^n ~ / ^ / ^ 
i=l j=i 

Xn+1 — OLn^n \ t-n+l 



J-j.n 



Otri. — i ' ' ' OL 



"''n— i+1 



n-j 



(33) 
(34) 



Proof. By Corollary [6] it is only necessary to determine a sequence \an\ of nonzero elements of J- 
sucli that for each n the quantity (|25|) is independent of uq for the following function 

/„(uo, . . . , -Ufc) = ai,„Uo + a2,n'Ul H h a^.n^fc + &n- 

For arbitrary uq, f i, • • • , Vfc S J-" and j = 0, 1, . . . , fc define Ci,n(ii0 5 'Wij • • • , Wj) as in Corollary [6l 
Then the expression (j25p is 



-On-uo + 6„ + ai,„no + a2,„Ci,n(i'o, ^i) H h afc,„Cfc,n(u.o,fi, ■ . . ,ffc) 

k 



bn + 



E 



Oi 



■J," 



t Ili=l "n- 



a.. 






The above quantity is independent of uq if and only if the coefficient of uq is zero for all n; i.e. 
if {a„} is a solution of the Riccati difference equation 



E 



-^j,n 



1 lli=l O^n-i 

which is Eq. (l32p . It follows that Eq. (l3ip has a linear form symmetry of type (I23p with eigensequence 
{«„} for each solution {a„} of the Riccati equation. For the corresponding SC factorization of 
()3ip . the cofactor equation is simply ()24p while the factor equation is obtained using the above 
calculations and Eq. (|20p of Theorem d] as follows 



tn+l — bn — 2_^ ^j,n 2_^ 
3=1 

k k 



^71— i+1 
1 1 [p=i '^n-p 



""-EE 

1=1 j=i 



^],n 



Oir. 



■ a. 



'■n— i+l' 



n~j 



This completes the proof. ■ 

Corollary [9] states that any solution of the Riccati equation ()32p gives a form symmetry and a 
SC factorization of (j3ip as specified above. The next example illustrates Corollary [9l 
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Example 10 Consider the second- order difference equation 

x„+i = (-l)"+ix„ + x„_i + &„ (35) 

where bn, xo,x^i are in a field T which we may take to be any one of the familiar fields Q, M or 
C. The associated Riccati equation of i35\) is 

a^ = {-l)^+^ + J—. (36) 

On-l 

Straightforward calculation shows that if uq ^ 0,-1 then 

oo + 1 1 

ai = , a2 = — r, as = -ao, 

ao ao + 1 

ao + 1 1 

a4 = , as = — -, ae = ao- 

ao ao + 1 

It follows that all solutions of the Riccati equation jg6]) with initial value outside the singularity set 
{0,-1} are eigensequences in J^ of period 6: 

ao + 1 1 ao + 1 1 

"0, , — r,-ao, , — r,ao, • 

ao ao + i ao ao + i 

The SC factorization of the linear equation \35^) is now obtained by Corollary O as 

a„-i 

The next result is concerned with the case of constant coefficients. The straightforward proof 
is omitted. 

Corollary 11 Let {6„} be a given sequence in a non-trivial field T and let {oi}, i = 1, . . . ,k be 
constants in T such that at 7^ 0. 

(a) The non-homogeneous linear equation of order k + \ 

Xn+l = aoXn + aiXn-l H h akXn~k + bn (37) 

has a linear form symmetry with eigensequence {an} for every solution {an} in J- of the following 
autonomous Riccati equation of order k 

an = ao-\ \ \ \ . (38) 

a„_i a„_ia„_2 a„_i • • • a„_fc 

(b) Every fixed point of l[38\) in J^ is a nonzero root of the characteristic polynomial of |37[ ), i.e., 

A*-'+i - aoA'^ - aiA^-i Ok-iX - Ok (39) 

and thus, an eigenvalue of the homogeneous part of I^Tj ) in T . As constant solutions of i38\) such 
eigenvalues are constant eigensequences of (37^ . 
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Example 12 Consider the autonomous second- order linear difference equation 

Xn+l =Xn + Xn-1 (40) 

Eg. 1^0^ has two real eigenvalues 

as roots of the characteristic polynomial A^ — A — 1 or equivalently, as fixed points of the Riccati 
equation 

an = l + -^ (41) 

(Xn-l 

in the field J- =R. Thus each of a+ and a_ is a constant eigensequence of I^Oj ) in M and the 
following SC factorization is obtained in M; 

_ 1 _ 

Xn+l = a+Xn + tn+1- 

Note that the SC factorization above has constant coefficients also. We note further that since a± 
are irrational the above SC factorization is not valid if T =Q the field of rational numbers. In 
fact, since the characteristic polynomial has no rational roots, it follows that there are no constant 
eigensequences for |^0[ ) in Q. However, Riccati equation ^4^\ ) ^■5 a rational equation and thus with 
a rational initial value oq the corresponding solution of |^i[ ) is a solution (non-constant) in Q. For 
instance, if oq = 1 then the corresponding solution of ( f^ip is an = (fn+i/v'n where {^Pn} is the 
Fibonacci sequence 1,1,2,3,5,8,. . . This rational eigensequence yields the following SC factorization 
of \4C^ that is valid in Q; 

_ ^n-l 
''n+l — ''Ji; 

^n 

_ V^n+1 
Xn+l — Xn \ iJn+1- 

^n 

We note that Ivca.n^^ipn+i/'-Pn = ck+ ii^ the above factorization; in this way the factorization over 
rationals is related to the earlier factorization over the reals. In a similar fashion, the equation 



has two complex eigenvalues 



Xn+l — Xn Xn—1 V^-^j 

l±iV3 



a± 
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that are roots of }? — A + 1. Thus, |^g[ j has no constant eigensequences in M but it does have 
non-constant real eigensequences since the Riccati equation 



"n = 1 



"n-l 



with the initial value ao = 2 has a solution 



2,1,-1.2,1,-1. 

'2 2 



of period three in R with a corresponding real SC factorization 

1 

/n contrast to the factorization of Eq. |^0[ j i/iere is no simple relationship between the factorization 
of (4W over the real numbers and its factorization with constant eigensequences over the complex 
numbers. 

Remark 13 Is it possible that a linear difference equation has no eigensequences, constant or 
otherwise in a given field J- because the associated Riccati equation has no solutions at all in Tl 

We know the answer to this question in some cases. If we have a linear equation (homogeneous 
or not) with constant coefficients in an algebraically closed field T (e.g., the field C of complex 
numbers) then T always contains constant eigensequences, namely, the roots of the characteristic 
polynomial Ii39\) . On the other hand, for the finite field Z3 = {0, 1, 2} with addition and multiplica- 
tion defined modulo 3, the linear equation ^0) has no eigensequences. This can be shown by testing 
each of the two possible nonzero initial values 1,2 in the Riccati equation ^1^ to verify that both 
lead to the singularity at 0: 

ao = 2^ai = l + - = l + 2 = 0, 

The answer to the question of existence of eigensequences in the general case is not known at 
this time; in fact, it is not known if a linear equation with real coefficients exists that has no real 
eigensequences. For "large" fields such as M or C it seems likely that the general linear equation 
/i31\) has an eigensequence in the field. 



The occurrence of Riccati difference equation in Corollary [9] may seem less surprising if we recall 
some basic facts from [5j. In particular, the homogeneous part of ()3ip is a homogeneous equation of 
degree one relative to the multiplicative group J-"\{0}. Therefore, it has an inversion form symmetry 
and the factor equation of its SC factorization is none other than the Riccati equation (I32p . Using 
this fact it is possible to restate Corollary [9] without explicit reference to the Riccati equation as 
follows. 
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Corollary 14 Assume that the homogeneous part of Eg. \31]) has a solution {i/n} in the field T 
such that yn 7^ for all n. Then {yn+i/yn} is an eigensequence of i31\} whose SC factorization is 
given by the pair of equations 13'5'j) and p'^ . 

Proof. It is given that {yn} satisfies the homogeneous part of (j3ip . i.e., 

Vn+l = aO,nyn + Ol.nJ/n-l + «2,nyn-2 + • • • + ak^nVn-k- 

Since yni^ ^ for all n, we may divide the above equation by y„ to obtain 

Vn+l , Un-l , Vn-2 , , Vn-k 

— OO.n + ai,n \- a2,n !"•••+ ak,n 

Un Un Vn Vn 

Vn-l , yn~2yn-l , , Vn-k Un-l 

— CiOn + ain i- a2,n !-••• + a^n • 

Vn Vn-l Vn Vn-k+l Vn 

Now defining a„ = yn+i/Vn for all n and substituting these terms in the last equation above 
yields the Riccati equation (i32l) . Thus {yn-^i/yn} is an eigensequence of (i3T]l in J^\{0}, as claimed. 
The SC factorization is obtained as in the proof of Corollary [9j ■ 

Corollary 15 In Eq. h31\) let {ai^n}, i = 1, . . . ,k and {6„} be sequences of real numbers with ai^n ^ 
for all i, n and Ok^n > for all n. Then ^31\) has an eigensequence {yn+i/yn} o,nd a SC factorization 
in M given by the pair of equations ^3^ and (34\ )- 

Proof. If we choose y-j = 1 for j = 0, . . . ,k then the corresponding solution {vn} of the homo- 
geneous part of (|3ip is a sequence of positive real numbers. Now an application of Corollary [T3] 
completes the proof. ■ 
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